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Abstract 

Using hadron dynamics we calculate nuclear transparencies for protons, 
knocked-out in high-Q^, semi-inclusive reactions. Predicted transparencies 
are, roughly half a standard deviation above the NE18 data. The latter 
contain the effects of binned proton missing momenta and mass, and of fi- 
nite detector acceptances. In order to test sensitivity we compare computed 
transparencies without restrictions and the same with maximal cuts for miss- 
ing momenta and the electron energy loss. We find hardly any variation, en- 
abling a meaningful comparison with data and predictions based on hadron 
dynamics. Should discrepancies persist in high-statistics data, the above may 
with greater confidence be attributed to exotic components in the description 
of the outgoing proton. 
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The following note concerns the transparency of a nuclear medium for a proton, knocked- 
out in a high-Q^, semi-inclusive (SI) reaction. Simplest is the electron-induced SI reaction 
A{e, e'p)X, where in addition to the scattered electron, one also measures the kinematics of 
a knocked-out proton. For it, the transparency T = T'^^ is the ratio of the experimental 
yield and some reference cross section. Desiring it to be a measure for the interaction of 
the exiting nucleon with the remaining nuclear core, the natural choice for reference cross 
section is the A{e, e'p)X yield where that interaction is absent, i.e. the (non-measurable) 
cross section in the Plane Wave Impulse Approximation (PWIA) Also in use is Z times 
the ep Mott cross section, but in that case the resulting ratio does not tend to 1 in the above 
limit. 

To the extent that the knocked-out particle is a hadron, conventional nuclear dynamics 
should suffice for a description of T. However, for sufficiently high energies the description 
of the knocked-out particle is predicted to require sub-nucleonic degrees of freedom, which 
interact only weakly with the nuclear medium 0. A feducial detection of the resulting color 
transparency (CT) from high-Q^ knock-out reactions then clearly requires 

i) A measurement of T under well-controlled conditions. 

ii) An accurate calculation of the non-exotic component of the transparency. 

Both requirements are hard to meet for SI (e, e'p) processes. Data are in general not 
genuinely semi-inclusive, but contain contaminations, e.g. cuts in observables and detector 
acceptances On the theoretical side the complexity of the many-body problem stands 

in the way of an accurate calculation of T. In the absence of deconvoluted data, calculations 
have sometimes modeled the above problematic features, relaxing to some extent the measure 
of semi-inclusivity. This has led to confusion and it is not always clear what is meant by 
'a standard Glauber calculation' for a SI transparency. For clarification we shall use the 
approach of Refs. with modifications expounded in 

Having defined the transparency T and the framework for computation, one may use the 
same to calculate related quantities, which experimental conditions may force one to study. 
Those become some measure for T, but are not identical to it. The spread in outcome of the 
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analysis then establishes bounds for the sensitivity of T for various, controlled deviations 
from semi-inclusivity. 

Consider the SI yield when the electron transfers momentum and energy (g, uj) to the 
target and the knocked-out nucleon has momenta p. Having factored out the (off-shell) 
Mott cross section from the yield, the focus is on the remaining SI response per nucleon |^ 

S''iq,uj,p) = l/AY,\Ton{q)\'S{u + A„ - ep) 

n 

^on(g) = {K\p\\Kp) (1) 

Here jFo„ is the inelastic charge form factor for a transition between the ground and an 
excited scattering state. The latter describes asymptotically a proton with mass M and 
total on-shell energy e{p) = {p^ + M^)^/^ and a core oi A — 1 nucleons in an excited state 
n, separated from the ground state by an energy Exploiting the high momentum of 
the exiting proton one may approximately factorize those into states for the core and the 
knocked-out proton scattered in the field of fixed scatterers 

Kp ^ *n"VJ (2) 

Substitution of (2) into (1) and subsequent application of closure leads to [^|| (s = — r^) 



Tfl if ~- 

S^\q,yo,p) = j5{yo-Pz) J dn J dse'P^ pi{ri,r[)R{q, r[,s. 
The energy loss uj in (0) has in ( pa] ) been replaced by a scaling variable 

yoJ V 



yo ~ yo 
yo = 



1 

2AV 



M 



q + ^2M{uo - {A)) + {u - {A)y 



(3a) 
(3b) 



(4) 



with (A), an average nucleon separation energy for the nuclear species considered. 

For any inclusive process the nuclear input required in ( P5| ) is the nuclear density matrix 
PAi diagonal in all coordinates except in the one of the struck particle ('!'). In contradis- 
tinction, the expression for 7 in (|3b|) which describes Final State Interactions (FSI) depends 



on the type of the SI process as well as on the approximation used in the evaluation. The 
induced high-energy rescattering of the knocked-out nucleon from the core is for all processes 
routinely described by Glauber theory p!0 |. 

One thus finds from (0), (0) for the 'unrestricted' SI process (i.e. without binning of 
data) § 

7(g, r, s) = (l + Tff{b, z)) (l + [rf z + s)]*) , 



(5) 

where for short-range NN interactions V one may approximately relate the elastic off-shell 
profile function F"-^-^ in (^) to its corresponding on-shell analog F. With Vq = e{q)/q 



Ff ^(6, z) ^ exp[-{z/v,) / d(V{b, ()] - 1 ^ exp[-iz/vq)e{-z) / dCV{b, ()] 



= e{-z)Tq{b) = e{-z)(exp[-{z/v,) j^jcv{bx)] - 1) 
For the on-shell profile F we use the standard parametrization 

Vqib) = e'^^^"^ - 1 ^ -(af72)(l - tTq)A{q, 6), 



(6) 



(7) 



with r, the ratio of real and imaginary parts of the forward A^A^ elastic scattering amplitude, 
and A{q,b) accounting for the range of Fg(6); zero-range corresponds to A{q,b) —>■ S'^{b). 
Inelastic contributions to elastic pp and pn scattering are accounted for by the imaginary 
part of the eikonal phase x (0) ^^"^ impact parameter representation related to 

the 'partial' inelastic cross section or^'*"'^^(6) 



Using (g), 
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p,inel 



(6)^1 



-21mxg(fe) 



Aiq,b)a 



inel 



A(g,6)^^exp[-(6g72)2] 



Air 



a 



inel 



„tot el 
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one shows that 7 in Eq. (FT) becomes 



7(g,r,s) = -A{b) 
{(3a) = {A-l) 



-z){a''"'^) + e{z)9{s- z) 



l + ir 



-cr 



tot 



-1 



(8) 



(9) 



where we defined averaged (weighted) pN cross sections. 

Eqs. (H) describe in principle FSI to all orders in 7, but for large momentum transfers 



it suffices to retain terms up to second order in 7. Those involve p„, n < 3 p]l|JT2 

rn-l 



i>2 



Cn{ri,rj;s) 



(2(^-1, s) ^ ^Jg{ri - r2)g{ri - ra + s,q) 
C„(ri, .., s) ^ [| (2{rj, r^, s^) (10) 

In all pn we factored out the non-diagonal single-particle density pi{ri,r[) which may ap- 
proximately be parametrized as |13[ {S = (ri + r'^)/2; s = ri — r'^) 



Pi 



:r^,r[)^piS) J dS' p,is, S') ^ p{r^)J:is) 



and where p(r) = pi{r,r) is the single nucleon density. Next one encounters in the ex- 
pression for P2 in (|TOD a non-diagonal pair-distribution function (2, approximated there by 



means of the diagonal pair-distribution function g |T2[. The non-diagonal n-body distribu- 
tion function (n in ( P^ is, in the independent pair approximation expressed as a product of 
pair analogs ^2- Using (|^) and ([TlD , Eqs. @) become 



S''{q,yo,p) = "^Sivo-p.) J dse'P^^{\s\) J dr,p{r,)R{q,r,, s,) 

= -S{yo - P.) I dse^P^i:{\s\)G{q. s^) (12a) 
q J 

G{q, s,) = [ dnp{r,)R{q, n, s,) (12b) 



5^™^(g,yo,p) = ^5(yo - P.) / dse^P^^{\s\) J dr,p{r,) = {m/q)6{yo - pMp) (12c) 
Treating the retained FSI terms in the first cumulant approximation, one finds for A ^ 1 

R = e^ = e^2+f^3+... (13a) 
^2iq, ri; s) = (A ~ 1) J dr2p{r2)l{q, ri - rs, s)C2(t*i - ^"2, s) (13b) 
^3{q,ri;s) = J y c?r-2c/r3p(r2)p(r3)7(g,ri - r2,s)7(g,ri - r3,s)C3(T'i,r-2,r3,s) (13c) 
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We shall later on return to the ternary collision terms above and discuss first the binary 



collision part in (|T^) P | 



■ytot Tyinel 



1X2 — ^2 



R^iq; ri,s,) ^ exp 



(fbA{b) / dzp{bi — b, zi — z)(2ib, Sz 



(14a) 



Ri''^\q; ri,s,)^ exp 



(A - 1)^(t;"^'^ J d'bA{b) 
/ dzp{bi-b,zi- z)C2{b,z,s, 

J —oo 



(14b) 

We note that the results Eq. (|12|)-(|14|) contain measurable quantities without reference 
to an underlying potential method. It seems then reasonable to dissociate the outcome from 
that model and to postulate its validity also in the high-Q^ regime. 

We start with Eqs. (|12aD, (|12cD for the 'unrestricted' nuclear transparency, without any 



binning or correction for detector acceptance 



Tiq,p)^ 



^^^'™(g,l/o,p) 



J dseP^i:{\s\) 



(15a) 



mip)]-' I dse^P^J:{\s\)Giq,s.) = [nip)]-' [ M^(p^,p, - p'JG(g;p'J (15b) 
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T is from (1) expected to be a function of q, the energy loss u or i/q, Eq. (^, and the 
momentum of the knocked-out proton momentum p. The fact that (|I3p does not appear 
to depend on is an artifact of closure. Without its application one is lead to ([L5b|) with 



non-canceling single-nucleon spectral functions instead of momentum distributions. 

Next we consider conditions, intermediate between semi and total inclusive scattering 
and consider first cuts on the missing momentum = P — Q- For fixed q those are the 
same as cuts on p. We thus integrate both SI yields or responses (^ in the ratio (|15a]) over 
the momentum p of the outgoing proton, with the result 

/ dpS'^iq, yo, p) J ds,e^yo^^J:{s,)G{q, s,) 



T^(g,l/o) 



JdpS''^i'P^^\q,yo,p) 
dp' 



J ds.e'yo'^Eis,) 



{2n) 



Fo{yo-p')G{q,p') 



(16a) 
(16b) 



and where we used the asymptotic hmit of the NR total inclusive response [^2 

dP 



Foivo) = lim —S'^\q,yo) 

q^oo q 



(17) 



After the remark following Eq. (|T3p one observes the reappearance of the energy loss through 
yo in T^, Eq. (|IBp. Notice that the response S^^'^^^^ in the ratio in ( |16aD , when inte- 
grated over the the proton momentum and energy loss is 1, which implies Zct^°** for the 
corresponding reference cross section. This does not contradict the original choice of PWIA 
yields as reference for the definition of T: it emerges if in (|15a|) one consistently applies the 
same cuts in actual and reference yields. 

Finally we consider the transparency when the cross sections in the ratio ([T6|) are in 
addition integrated over the energy loss of the electron. From Eqs. (01) and (M)) one readily 



finds 111 



r^^(g) ^ Giq, 0) = j dr,pir^)R{q, r,, 0) 

dr,p{r^)exp[-{A - J d^bA{b) dzp{h, - b, - z)g{b, z)] 



U-l){g){a-')]-' 



db 



1 _ exp( - (A - l){g){al"^'){h{b) + h{b) + ..} 



(18a) 
(18b) 
(18c) 



Eq. ( |18c| ) holds only in the 0-range interaction limit and t2{b) = p(6, z') is the standard 
thickness function resulting from binary FSI. ^3(6) there is due to ternary contributions, to 
be discussed shortly. 

Although the actual derivation above leads to o"*"^', one occasionally meets expressions 
for T with o"*""^' 



a 



tot 



where in particular the binary part 



dr^p{r^)exp[{-iA-l){af') dzp{h,, z, - z)g{Q, z)] 



U-l){g){af')]-' 



d'b 



l_exp -(A-l){g){a'"')t2ib) 



(19a) 
(19b) 



is reminiscent of the transparency for an elastically scattered proton. The above has for 
instance been used by Benhar et al, who extended a previous description of the totally 
inclusive (e, e') process ||2^ and reached T^^, Eq. ( |19a| ) [ |21|| . The same holds for Frenkel, 
Frati and Walet and the criticism voiced in and in particular in pertains to it 



as well. Finally we mention work by the Rome-Perugia group, who use a spectral function 
beyond the mean-field approximation, but do not seem to go beyond the PWIA p3[ . 

We now discuss Eq. ( |18c| ) in some detail. Its binary collision version (with ternary and 
higher order corrections ts + ...—>■ 0) resulted when Kohama etal [^,0 applied unnecessarily 
strong approximations in an evaluation of the 'unrestricted' transparency T, Eq. (P^). 

Alternatively, desiring to simplify matters, Nikolaev etal just replaced T T^^, ob- 



tained from cross sections integrated over missing momentum and energy- loss [|l^]. Those 
authors considered binary ('hole') as well as ternary ('spectator') contributions with the 
scattering operator ( p.8a| ), appropriate to T^^ , i.e. (§) with s — > 0: The integration over the 
energy loss eliminates the retardation in the propagation of the density disturbance. As a 
consequence static (diagonal) quantities replace everywhere non-diagonal ones, e.g. (2 9, 
etc. 

Instead of g Nikolaev et al prefer in ([T0|) the use of the pair-correlation function C = 1 —g, 
thus 

^3(1, 2, 3) = (3(1, 2, 3, 0) ^ g{l, 2)g{l, 3)g{2, 3) (20a) 
«C(1,2)C(1,3)-C(2,3) (20b) 

where for brevity we wrote g{12) = g{ri,r2), etc. Making the additional approximations 
g = —C and the second part of ( pOb|) in Eqs. ( |13b|) , respectively (|13c| ) for s = 0, those 
become 

^r'^^lg, r,) ^ (A - 1) I dr2p(r2)7^^(g, r, - r2)C{r,, r^) (21a) 
^r^'^iq, ri) ^ ^ / / dr2dr,p{r2)p{rsh''''{q, r, - r^h'^'^iq, r, - r,)C{r2, r,) (21b) 

Next one shows that the product CC in ( |20b| ) equals (1/2)^^2^-^'^^ and similar approxi- 
mations to higher order densities produce terms of the same order in the expansion of the 
first cumulant of ( |21aD . Nikolaev etal thus reach 

R"" = exp(n^ {[C]) + + 

Q^[C])=Q2[9]) (22) 
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In view of the replacement g —C in the actual expression Eq. (|13b|) for ri2, that approxi 



mation misses contributions to Vt^ linear in 7. Likewise, there are ternary FSI contributions 
quadratic in 7 which are not contained in either 1/2 (fi^)^ or in ( pib| ). 

Comparison of numerical results reflects the above shortcomings. Nikolaev etal report 
that the first cumulant of binary as well as ternary FSI contributions Eqs. (^Tj) are + (5- 
6)%, respectively -(2-3)% of the PWIA result. However, the former fraction, computed from 
( |13b|) with C —g is much larger, namely 15-20 % ||T9|,[7[|. Even assuming (??) to be a 
measure for ternary FSI, those cut the binary part by only ^ 15% and not by ~50%. Clearly, 
in view of the approximations involved, it seems safe to disregard ternary FSI altogether. 

We now present some numerical results for the transparencies of a number of nuclear 
species under NE18 kinematics 0,^, i.e. for Q'^=1.0, 3.0, 5.0 and 6.7 GeV^, p ~ \p\q. For 
the evaluation of the predictions Eqs. ([12|)-(p!8D one needs nucleon densities p(r) and pair 
distribution functions g (both taken to be the same for protons and neutrons) and dynamical 
input. Measured or interpolated values for representative A^A^ scattering parameters are 
given in Table I. A glance at the entries suffices for one prediction: In the region of the 
NE18 experiment transparencies should, for any given nuclear species decrease as function 
of and reach a plateau beyond ~ 1.5 GeV^. Any significant deviation from that 
prediction cannot be accounted for by standard hadron physics. 

In Figs. la,b,c we entered for C, Fe and Au as function of the undeconvoluted NE18 
data 1^ and those are compared with the following predictions (none contain corrections for 
finite solid-angle detector acceptance). 



i) Unrestricted transparencies T, Eq. (|T^), caused by binary collisions, which in turn 
are generated by finite, as well as by zero-range interactions. 

ii) Results for (p!6D, when cross sections have been integrated over missing momentum. 
Compared to T, changes are < 2 %, independent of the nuclear species or Q^, and are 
insensitive to any reasonable choice for the average separation energy. 

iii) T^^ , Eq. ( |18a| ), obtained from cross sections, in addition integrated over the electron 
energy loss. Relative changes now amount to (4-6)%. We checked that in all cases only tiny 
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changes occur if one discards the non-neghgible r (Table I). Out of all parameters studied, 
only the A^A^ interaction range affects predictions. 

iv) Sizably smaller transparencies T^^, Eq. ( |19b| ), result from Eq. 
cross sections replace the smaller inelastic ones. 

We conclude: 

1) For all targets and all Q^, the most complete calculations using binary FSI produce 
transparencies, roughly half a standard deviation higher than the data. Lacking data with 
better statistics, it seems premature to ascribe the current discrepancies to color trans- 
parency, whose onset is anyhow not expected to occur for as low as ~ 1-1.5 GeV^. 
Ternary FSI contributions decrease the above results by only a few percent. 

2) There are only modest changes in T due to acceptable ranges in input parameters. 
Likewise T is only marginally sensitive to 'maximal' cuts and it seems reasonable to assume 
the same for actual, more restricted cuts. 

3) Having demonstrated the stability of calculations based on conventional hadron dy- 
namics, firm conclusions can be drawn regarding the existence of exotic parts in T, should 
discrepancies with high-statisitic data persist. Presently one can only say that no conven- 
tional hadron theory can predict the rise in the undeconvoluted NE18 data towards the 
largest measured points for all targets, or the relatively flat carbon data for the lowest 



19a), where total 
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Figure Captions. 

Fig. la. Transparencies of C tlie passage of a proton, knocked out in a semi-inclusive 
A{e,e'p)X reaction for NE18 kinematics as function of Q^. Data are from O'Neill etal ^ 
and differ slightly from previously published data by Makins etal The drawn, dotted 
and long dashed lines correspond to respectively T, T^, T^^ , Eqs. (0), (|T§) and ( p.8c[ ). 
The dot-dashed line gives the 0-range interaction result for T. Short dashes are for T^^, 
Eq. (|l9iD . 

Fig. lb. Same as Fig. la for Fe. Data are from O'Neill etal 
Fig. Ic. Same as Fig. lb for Au. 



Table I. Partly interpolated pp and pn elastic scattering parameters 24 



(GeV) (GeV^) 


pp pp 

(mb) (mb) 


^tot ^inel 
pn pn 

(mb) (mb) 


'^pp '^pn 


{Qq)pP {Qo)pn 

(GeV) (GeV) 


1.20 1.04 


36.6 12.6 


36 8 


0.30 -0.25 


0.84 0.84 


2.40 3.0 


46 26 


42 28 


-0.35 -0.40 


0.59 0.59 


3.54 5.0 


42.4 27.1 


44.0 27.1 


-0.37-0.45 


0.56 0.56 


4.48 6.7 


42.4 27.5 


44.0 27.5 


-0.37-0.45 


0.56 0.56 
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